Modelling of compound nucleus formation in fusion of heavy nuclei 
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A new model that includes the time-dependent dynamics of the single-particle (s.p.) motion in 
conjunction with the macroscopic evolution of the system is proposed for describing the compound 
nucleus (CN) formation in fusion of heavy nuclei. The diabaticity initially keeps the entrance 
system around its contact configuration, but the gradual transition from the diabatic to the adiabatic 
potential energy surface (PES) leads to fusion or quasifission. Direct measurements of the probability 
for CN formation are crucial to discriminate between the current models. 
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The understanding of new experimental results on fu- 
sion of heavy nuclei and the formation of superheavy el- 
ements (SHE) require modelling, not only of the initial 
capture process and the final CN de-excitation process, 
but also of the intermediate stage of evolution of the com- 
bined system from the contact configuration into the CN. 
The competition between fusion and quasifission (resep- 
aration before CN formation) can inhibit fusion by many 
orders of magnitude, e.g. Understanding this inhi- 
bition may be the key to forming more SHE. Nowadays, 
there is no consensus 2] for the mechanism of the CN 
formation in fusion of heavy nuclei near the Coulomb 
barrier. Depending on the main coordinate for fusion, 
two sorts of models can be distinguished. In the first 
type mmillllilillillll, the fusion occurs along 
the radial coordinate using adiabatic PES obtained either 
with the liquid drop model or Strutinsky's macroscopic- 
microscopic method. The competition between fusion 
and quasifission, which depends on the fluctuactions Ra| , 
has only been included in recent models, e.g. j9l[Tolll]|. 
However, the experimental data are not always explained. 
In the second type 13] (dinuclear system (DNS) model), 
the fusion happens in the mass asymmetry coordinate 
r] = {Ai — A2)/ {Ai + A2) where Ai and A2 are the mass 
numbers of the nuclei. The DNS nuclei remain at the 
touching configuration and exchange nucleons until ei- 
ther all nucleons have been transferred from the lighter 
to the heavier fragment (complete fusion), or the DNS de- 
cays before the CN formation (quasifission). The model 
assumes a sudden (double-folding in frozen density ap- 
proximation) PES in the radial coordinate, while the PES 
behaves adiabatically along the fusion path in the rj co- 
ordinate. Although the model has been used to explain 
many experimental evaporation residues (ER) cross sec- 
tions, its theoretical foundation is not clear enough yet. 

In this paper a new model is proposed for the CN for- 
mation, which is based on the following general ideas 
(which are well-established but have up to now not been 
used in combination in any of the current models of fu- 
sion): (i) Once the two nuclei are at the contact point, 
the system moves in a multidimensional space of collec- 



tive coordinates q (shape parameters), (ii) this motion is 
governed by the master equation, and (iii) the nature of 
the s.p. motion is time-dependent, it is initially diabatic 
and then approaches the adiabatic limit due to residual 
two-body collisions. Consequently, the system moves in 
a time-dependent PES, V{q, t), which is initially diabatic 
and gradually becomes adiabatic. In the diabatic limit 
[T^ (elastic nuclear matter), the nucleons do not occupy 
the lowest free s.p. energy levels as in the adiabatic case 
(plastic nuclear matter), but keep their quantum num- 
bers and remain in the diabatic states during a collective 
motion of the system. This approach is realistic in the 
initial stage of collisions near the Coulomb barrier where 
the total excitation energy per nucleon E* > 0.03 MeV 
ITgII . During the transition from the diabatic to the 
adiabatic limit, the nuclear matter is elastoplastic like 
glycerine. In contrast to the current models for CN for- 
mation, the present approach explicitly includes for the 
first time the time-dependent dynamics of the s.p. mo- 
tion in conjunction with the macroscopic evolution of the 
system into the CN. In Rcfs. 0, Q, ^ JOiJ the ideas (i) 
and (ii) were applied, but the authors exclusively used 
adiabatic PES. In a recent paper H^, only idea (iii) was 
used to calculate the dynamical potential for the radial 
motion of the combined system, while in the past only 
the diabatic limit of the s.p. motion was used to study 
the initial capture process 0, 0|. The ideas (i)-(iii) 
have solely been applied in combination to describe deep- 
inelastic reactions [l^- The scenario presented here for 
the CN formation shows that following contact the di- 
abaticity forms a valley in the PES where the system 
remains trapped around its touching configuration, but 
the gradual transition from the diabatic to the adiabatic 
PES allows the system to evolve in shapes leading to fu- 
sion or quasifission. The timescale for the decay to the 
adiabatic PES is crucial in determining the timescale of 
the transition to the compact fused system. The calcu- 
lations are based on the master equation |20l and the 
diabatic two-center shell model (TCS M) [211 developed 
using the asymmetric TCSM (ATCSM) p^The critical 
ingredients of the model are (i) temperature, (ii) diabatic 
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PES, (iii) adiabatic PES, (iv) transition from diabatic to 
adiabatic PES, and (v) shape parameters. The evalua- 
tion of each will be described, then the results from the 
model will be presented. 

The time-dependent population probabilities p for the 
different configurations q (shapes) of the system are so- 
lutions of the master equation 



the adiabatic PES which is calculated using Strutinsky's 
macroscopic-microscopic method |22| . The nuclear part 
of the macroscopic energy is obtained with the Yukawa- 
plus-exponential method '25]. The diabaticity destroys 
the Fermi-distribution of the occupations, but the resid- 
ual two-body collisions gradually recover it. This process 
is locally described by the relaxation equation j26j 



p(q,i) =^[A(q,q',t)p(q',t) - A(q',q,i)p(q,i)], (1) 



where the macroscopic transition probabilities according 
to Ref. [13 are A(q,q',i) = Koea;p[F(q',<)/2r(q',t) - 
F(q, i)/2T(q, t)], justified by the assumption that the 
level density of the system determines the transition. The 
strengh constant kq characterizes the global time scale 



and has a realistic value of 



10^ 



The sum 



in Eq. ^ is extended only to the nearest configura- 
tions q (the collective coordinate space is discretized). 
It is assumed that the system is initially at the con- 
tact configuration qp where the s.p. occupation numbers 
obey a Fermi-distribution n^(qQ,To) for a temperature 
Tq. Configurations other than qg are not populated at 
this time, and hence the initial condition for Eq. iQJ is 



p(q,o) 



The temperature Tq is related to the 



excitation of the system immediately after the capture 
process. This temperature can be estimated either as 
To w yJ[E(..m. — ^(qoiO)]/a, where Ec.m. is the total in- 
cident energy in the center of mass frame and a = A/ 12 
MeV~^ {A is the total mass number of the system), if the 
initial radial kinetic energy is dissipated when the nuclei 
reach the contact point, or using a frictional model, e.g. 
[To| . for the capture process. The local (at fixed q) tem- 
perature T = y^Eexc/o, is defined by means of the local 
excitation energy 



dAVdiab{q,t 
dt' 



■)pici,t)dt , (2) 



which results from the decay of the diabatic part AVdiab 
of the potential V. AVdiab represents an energetic hin- 
drance for the initial system to reach a configuration q, 
if the nucleons follow diabatic levels during this process 
(elastic response) . The local excitation of the system is 
caused by the loss of its elasticity 0|. AVdiab is calcu- 
lated as 



AVd^abiq, t) 



(q)K(q,i)-nf(q,r)], (3) 



where 

gdmfc are the diabatic levels with occupations 
Ua (q, t) and a denotes the quantum numbers of these 
states. The diabatic levels ef*"'' and their wave func- 
tions (r) are obtained with the maximum symme- 
try method The dynamical PES is defined as 
V{q,t) = Vadiab{<l,T) + AVd„fc(q, i), where Vadiab is 



-i(q,i)K(q,t)-nf(q,T)], (4) 



where r is an average relaxation time (in order to con- 
serve the number of particles). The initial occupations 
f^a(q, 0) are the diabatic ones obtained from n^(qQ,To). 
T is defined as 



EJ»„(q,0^<(q,T)]r„(q,r) 



(5) 



where F^ are the widths of the s.p. levels. The fac- 
tor Ncou is the average number of two-body collisions 
per nucleon to establish the equilibrium occupations n^. 
The value NcoU = 3 [^j will be used. The expression 
© follows the idea that the relaxation process becomes 
slower when the occupations approach n^. If the equilib- 
rium was reached, the relaxation time would be infinite, 
i.e. the occupations would remain the same. The widths 
r„ are obtained with the parametrization given in Ref. 
[27|. Since the diabatic s.p. excitations occur around the 



Fermi surface, the values Fn 



0.061 MeV"^ for half 



saturation density and c = 20 MeV will be used |23| . 

The collective coordinates q are the shape parameters 
of the ATCSM (i) the elongation A = l/2Ro, which 
measures the length / of the system in units of the diam- 
eter 2i?o of the spherical CN and describes the relative 
motion, (ii) the deformation (3i = ai/bi of the fragments, 
defined by the ratio of their semiaxes, (iii) the neck coor- 
dinate e = Eq /E , defined by the ratio of the actual bar- 
rier height Eq to the barrier height E of the two-center 
oscillator, and (iv) the volume asymmetry of the nuclear 
shapes (equipotential shapes) ^ = {Vi — V2) / {Vi + V2) , 
where Vi and V2 are the volumes of the left and right 
regions divided by a plane at the necks between the frag- 
ments. The collective coordinate space is divided into 
three regions: (i) compact shapes around the spherical 
shape (fusion region), (ii) elongated shapes outside the 
initial contact configuration and beyond the Coulomb 
barrier (quasifission region), and (iii) intermediate shapes 
which could lead to fusion or quasifission. In the fusion 
region the physical mass asymmetry rj (defined like ^ but 
in terms of the masses which are calculated using the mi- 
croscopic density distribution p{r) ~ X]q "■a|<^?"''('')P ) 
reaches a minimal plateau, i.e. a maximal number of nu- 
cleons move (their wave-functions spread) in the whole 
volume of the system due to the decrease of the bar- 
rier Eq between the fragments. The fusion {Pcn) and 
quasifission (Pqf) probabilities are defined as the sum of 
the population probabilities p in the fusion and quasifis- 
sion regions, respectively. In addition the CN excitation 
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FIG. 1: Dynamical PES for ""Pd+^^Pd; (upper part) at 
t=0 s (diabatic PES), the square denotes the contact config- 
uration; (middle part) at = 4 x 10~^° s; and (lower part) 
the adiabatic PES. See text for further details. 



energy Ecn is defined as the average of the excitation 
energies Q of the different shapes in the fusion region. 

In the following the fragments are considered as spher- 
ical {(3i—l) and the neck coordinate is fixed at e — 0.75. 
With this value of e the Coulomb barrier of the diabatic 
potential for an initial system with Tq = MeV is close 
to the barrier of the double- folding potential For 
computational reasons the calculations are done using 
the coordinates A and ^ that are the relevant ones in the 
current models for CN formation. The non-linear set of 
equations (|1I5|I are solved by sucessive iterations using a 
small time step Ai — 10^^'^ s. The master equation ^ is 
solved on a grid (1 ^ A 1.8,|C| ^ 0.7 where AA = 0.02 
and A^ = 0.1) with appropriate boundary conditions 
to follow the continuous split of the initial population 
probabilities p(q, 0) into fusion and quasifission. Values 
of 1^1 > 0.7 are not included in the calculation because 
the TCSM used is not appropriate for large ^. The 
fusion and quasifission processes are determined at the 
timcscale when PcN{to) + PQF{to) — 1- The model 
will be applied to some (near) symmetric central {I = 0) 
collisions. In these calculations the excitation energy of 
the initial system at the contact configuration is 40 MeV. 

Fig. 1 shows the dynamical PES for "°Pd-H"'^Pd (i) 
at the initial moment (diabatic PES, upper part) when 
the nuclei are at the contact configuration (square), (ii) 
at to = 4 X 10~^° s (middle part), and (iii) the adia- 
batic PES (lower part). The PES is normalized with the 
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FIG. 2: Distribution of the population probabilities p(A, ^, to) 
for "°Pd+"°Pd on the PES of Fig. 1 (middle part). The 
fusion region is A ^ 1.3. See text for further details. 



TABLE I: PcN, to (10"^° s) and Ecn (MeV) for some (near) 
symmetric central collisions. The Pqn values are compared 
to direct experimental P^^' [2^ . 



Reaction 


PcN 




to 


Ecn 


pcxp . 
^CN 


'JOZr+y"Zr^^«"Hg 


3x10" 


1 


5 


40 




iooMo+iooMo^2oop^ 


1.4x10 


-1 


4.5 


40 




110p^_^110pj^220u 


1.7x10 


-2 


4 


40 




^°°Mo+"°Pd->2iORa 


7.5x10 


-2 


4.3 


40 




««Zr+i24Sn^220^j^ 


4.5x10 


-2 


4 


40 


5x10"^ 



macroscopic energy of the spherical CN. In Fig. 2, the 
distribution of the population probabilities p(A,^,<o) on 
the PES of Fig. 1 (middle part) is presented. From these 
figures, it is observed that the diabaticity initially forms 
a valley confining the entrance system around its touch- 
ing configuration as in the DNS model 13]. Nevertheless, 
the gradual transition to the adiabatic PES leads to fu- 
sion mainly in the elongation A (relative distance) or the 
population of neighbouring configurations ending in the 
quasifission channel. In contrast to the DNS model, the 
diffusion along the ^ coordinate, while the nuclei remain 
at the contact point (A — 1.5 — 1.6), is strongly sup- 
pressed by a large diabatic hindrance. Pcn and Pqf are 
determined before the dynamical PES reaches the adia- 
batic one (comparing middle and lower parts in Fig. 1). 
The competition between fusion and quasifission is regu- 
lated by the dynamical PES. The probability of configu- 
rations which stay close to that of the entrance channel 
and reseparate is very large, as can be seen from Fig. 2. 
This causes the fusion hindrance for ^^"Pd-l-^^^Pd which 
is consistent with the experimental conclusions drawn in 
Ref. P]. The final population probabilities in the fusion 
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region (A ^ 1.3 in Fig. 2) are very small, and as a result of 
the expression Ecn is practically the same as the ex- 
citation energy of the initial system at the contact point, 
i.e. 40 MeV. The scenario seen so far is the same for other 
systems studied. Table I shows Pcn, to and Ecn calcu- 
lated for some (near) symmetric reactions. The values 
of the model parameters are the same for all the reac- 
tions. The decrease of the Pcn values from ^°Zr-|-^"Zr 
to ^^"Pd-l-^^^Pd is observed along with an increase of the 
quasifission probability beyond 90%. This is because (i) 
the repulsive character of the diabatic PES increases with 
increasing mass number A of the combined system and 
(ii) the valley of the diabatic PES becomes more shallow 
with increasing A. Consequently the timescale to of the 
fusion and quasifission processes decreases with increas- 
ing A. The Pcn for ^^Zr-|--'^^'*Sn agrees very well with 
the recent direct experimental measurement plj . The 
P^j^' includes the contribution of higher partial waves I, 
i.e. I > 0, although the main contribution results from 
low I for heavy systems at this energy. The Pcn obtained 
with the DNS model for this reaction ^-Iso agrees well 
with that experiment, while Pcn for, e.g. ^^"Pd-f ^"'^^Pd, 
is about one order of magnitude smaller than the present 
Pcn (see Table I). The fluctuation-dissipation model in 
Ref. Q predicts Pcn values that are much larger than 
the present ones, i.e. about one order of magnitude larger 
for ""Pd+"°Pd and i°''Mo-H"'^Pd. Direct experimental 
measurements of Pcn hke in Ref. are crucial to dis- 
criminate between the current models for CN formation. 
The dependence of Pcn on the parameter F^^^ is strong 
for the heaviest systems studied, i.e. changing the latter 
from 0.061 MeV"! to 0.03 MeY'^ causes the Pcn de- 
crease by about one order of magnitude due to the faster 
transition to the adiabatic PES, whereas the effect on the 



timescale is rather weak. The timescale to is mainly 
determined by the quasifission process that occurs near 
the contact configuration of the initial system, where the 
diabatic effects are small. The dependence of Pcn on 
the excitation energy of the initial system at the touch- 
ing point is rather weak, i.e. a saturation of the Pcn 
practically occurs from ~ 20 MeV upwards. 

In summary, a new realistic model for the CN for- 
mation in fusion of massive nuclei has been developed. 
It incorporates for the first time important physical ef- 
fects which critically affect the evolution of the fusing 
system. The diabaticity initially keeps the entrance sys- 
tem around its touching point, but the gradual transi- 
tion from the diabatic to adiabatic PES leads to fusion 
(mainly in the relative distance) or quasifission. The dy- 
namical PES regulates the competition between fusion 
and quasifission. The probabilities for CN formation in 
some (near) symmetric central collisions have been ob- 
tained and found to agree very well with the recent di- 
rect experimental determination for ^^Zr-|-^^^Sn. Direct 
measurements of the Pcn like in Ref. along with 

distributions of the quasifission fragments are crucial to 
discriminate between the current models for CN forma- 
tion. To calculate ER cross sections the present approach 
should be combined with other models that describe the 
initial capture stage and the survival of the CN against 
fission. A TCSM more appropriate for reactions with 
large mass-asymmetry is being developed which should 
be useful in predicting production cross sections of SHE. 
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